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We present a study of the results obtained combining LO partonic matrix elements with 
different orders of partons distributions. These are compared to the best prediction 
using NLO for both matrix elements and parton distributions. The aim is to determine 
which parton distributions are most appropriate to use in those cases where only LO 
matrix elements are available, e.g. as in many Monte Carlo generators. Both LO and 
NLO parton distributions have faults so a modified optimal LO set is suggested. 

The combination of the order of the parton distribution function (pdf) and the accom- 
panying matrix element is an important issue pQ . It has long been known that LO pdfs in 
some regions of x and Q 2 are qualitatively different to NLO (and NNLO) pdfs (see [2]) due 
to important missing corrections in splitting functions or coefficient functions for structure 
functions which are used to determine the 
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pdfs. Nevertheless, LO pdfs are usually 
thought to be the best choice for use with 
LO matrix elements, such as those avail- 
able in many LO Monte Carlo programs, 
though all such results should be treated with 
care. Recently it has been suggested that 
NLO pdfs may be more appropriate [3] , since 
NLO cross-section corrections are often small. 
There has already been an investigation of the 
use of NLO pdfs for the underlying event [4] . 
There is a big difference in the results when 
using CTEQ6L and CTEQ6.1M pdfs [5] due 
to the changes in the gluon, though agreement 
can be reached by significant retuning. This 
will affect predictions for other quantities. 

In this article we address the differences 
in predictions obtained for a variety of phys- 
ical quantities combining different pdfs with 
LO matrix elements. In each case NLO pdfs 
combined with NLO matrix elements repre- 
sent the best prediction - the truth. We inter- 
pret the features of the results and investigate Figure 1: Drell-Yan distribution at the LHC. 
how a best set of pdfs for use with LO matrix elements may be obtained. 

First, let us recall how LO pdfs tend to differ from NLO pdfs. The most marked dif- 
ferences are for light quarks at high x and the gluon distribution at low x. The coefficient 
functions for structure functions have ln(I — x) enhancements at higher pcrturbative order. 
This means the high-x quarks are smaller as the order increases. The quark-gluon splitting 
function P qg develops a small- a; divergence at NLO (with further ln(l/x) enhancements at 
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higher orders), so the small x gluon needs to be bigger at LO in order to fit structure func- 
tion evolution. Indeed, for Q 2 ~ 1 — 2 GeV 2 the NLO gluon is flat or valence-like at small 
x while at LO it grows much more quickly. These features show up in cross-sections. 

Let us start with the simple example of Drell-Yan production of vector-bosons from 
quark-antiquark annihilation as a function of rapidity at fixed invariant mass = 80GcV, 
i.e. appropriate for W, Z production, at the LHC. The NLO correction [5] for Drell-Yan 
production is quite significant, being positive for all rapidity and roughly 12% in this case. 
The absolute predictions and the ratios of LO matrix element results to the truth are show 
in Fig. [1] We see that indeed we are nearer to the truth with LO matrix elements and NLO 
pdfs [7] than LO matrix elements and LO pdfs [8]. Both LO ME results are too small, but 
NLO pdfs are closer and a better shape. LO pdfs and the LO matrix elements have the 
wrong shape being low at central rapidity but increasing at high rapidity where the high- a: 
quark is enhanced at LO. Hence, NLO pdfs seem more appropriate here. 

However, there is a small x counter- 
example. Consider the production of charm 
in low Q 2 DIS at HERA energies. The NLO 
matrix elements [5] contain a divergence at 
small x which is not present at LO. Using 
NLO pdfs the LO matrix element result is 
well below the truth at low scales and the 
shape is totally wrong, as seen in Fig. [2] The 
LO gluon is very large at small x since it has 
been extracted with missing P qg enhancement 
at small x, and this compensates the miss- 
ing small x divergence in the matrix clement 
- LO pdfs and LO matrix elements arc more 
sensible with a compensation between failings 
in both. Using LO pdfs and NLO matrix ele- 
ments the result is extremely large since there 
is a double counting of the small x divergence. 

From these two examples we can conclude 
that sometimes NLO pdfs are better to use 
if only the LO matrix elements are known, 
and we can get significant problems with the 
size and shape if LO pdfs are used. However, 
we can be completely wrong, particularly at 
small x, if we use NLO pdfs due to zero-counting of small- a; terms. Can we find some optimal 
definition of pdfs which have most desirable features? In order to make progress we need 
to better understand the difference between LO and NLO pdfs. The missing higher order 
terms in ln(l — x), l/x and ln(l/a;) in coefficient functions and/or evolution leads to pdfs 
at LO which are bigger at x — > 1 and at x — » in order to compensate. However, from the 
momentum sum rule there are then not enough partons to go around, and enhancements in 
some regions lead to depletion in other regions, particularly quarks at 0.001 < x < 0.1. This 
leads to a bad global fit at LO [8] - partially compensated by the LO extraction of as(M§ ) 
being ~ 0.130 to help speed evolution, and explains the underestimate of the Drell-Yan 
production at the LHC at more central rapidities. 

This obvious source of problems has lead to a suggestion [lOj that relaxing the momentum 
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sum rule could make LO pdfs rather more like NLO pdfs in the regions where they are 
normally too small. The resulting pdfs would still be bigger than NLO where necessary, i.e. 
the high-x quarks and low- a; gluon, but would not be depleted elsewhere. It is also useful to 
use the NLO definition of as. Because of quicker running at NLO, the LO and NLO couplings 
with the same value of as(M§) are very different at lower scales where DIS data used in 
global fits exists. Near Q 2 = 1 GcV 2 the NLO coupling with as(M'z) = 0.120 is similar to the 
LO coupling with as(M§) = 0.130. Hence, the use of the NLO coupling helps alleviate the 
discrepancy between pdfs at different orders. Indeed, the NLO coupling already used in some 
CTEQ LO pdfs [5] and in LO Monte Carlo generators. Relaxing momentum conservation 
for input pdfs and using the NLO definition of as does dramatically improve the quality of 
the LO global fit, X 2 = 3066/2235 for the standard LO fit becoming X 2 = 2691/2235, with 
a big improvement in the comparison to HERA data. The momentum carried by the input 
pdfs goes up to 113%. Using the NLO definition the value of as(M§) = 0.121. 

The modified pdfs, which we denote by LO*, become 
much more similar to NLO pdfs, in particular at small x 
the LO* quark distributions evolve as quickly as at NLO 
and arc similar for for x ~ 0.001 — 0.01. Similarly g(x, Q 2 ) 
is significantly bigger at LO* than at LO, and much bigger 
than NLO at small x. This will help when used with LO 
matrix elements for gluon-gluon initiated processes (e.g. 
Higgs production) where i-sT-factors are often much greater Table 1: a(bb) totals, 

than unity. We now look at the LO* pdfs in our first two examples, see Figs. [T]and[2] For 
Drell-Yan production at the LHC the LO* pdfs lead to shape of comparable quality as the 
NLO pdfs and the normalization is better. For the charm structure function comparing all 
possibilities the LO* pdfs and LO matrix elements are indeed closest to the truth. 

There is a similar conclusion 
for hadro-production of b quarks at 
LHC which probes low x ~ 0.001. 
We consider a(bb) where the ini- 
tial b quark has pt > 10 GeV 
and 1 77 1 < 5 using a NLO event 
generator [TT] and LO calculations 
[T2]. The total cross-sections are 
shown in Table. 1. The NLO 
pdfs and LO matrix element are 
clearly worst. We also illustrate 
final state px and rapidity distri- 
butions in Fig. [3] with pr > 18 
GeV for the b quark after shower- 
ing. Again the best absolute pre- 
dictions with LO matrix elements 
uses LO* pdfs and the worst NLO 
pdfs. However, in this case there is always a problem with the shape as function of pt- The 
NLO matrix element has a large positive effect at high px and very high rj. It is impossible 
for any parton shape to account for all NLO corrections. 

We also look at very high-Er jet production [13] at the LHC in Fig. [4] Ignoring the lowest 
Et, where hadronization and underlying event and possibly small x physics are an issue, 
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the LO and LO* pdfs with the LO matrix elements are of comparable quality to NLO pdfs. 
LO and NLO pdfs when combined with LO 
matrix elements produce results which devi- 
ate in opposite directions at high Et- Also 
shown is the prediction from LO** pdfs, 
which have an additional constraint for a 
harder high- a: gluon (this being small at 
LO even with momentum conservation re- 
laxed). This changes the results little com- 
pared to the use of the LO* pdfs. 

Hence, we can conclude that a fixed pre- 
scription of either LO or NLO pdfs with LO 
matrix elements or LO Monte Carlo gener- 
ators will lead each to incorrect results in 
some cases. To try to improve this situation 
we have suggested an optimal set of pdfs 
for LO calculations, the LO* pdfs, which 
are essentially LO but with various modifi- 
cations to make their features more NLO- 
like. These seem to work reasonably well 
and happen to achieve some of the features 
obtained by modifying pdfs in a process de- 
pendent fashion for use in Monte Carlo gen- 
erators/resummations discussed in e.g. |14j 
NLO matrix element corrections qualitatively change the features of the predictions for a 
physical process, regardless of how careful one is with pdfs, since new types of partonic 
process open up. This must always be borne in mind, and accounted for if possible. 
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Figure 4: The dijet cross-section at the LHC. 
More study is underway. However, sometimes 
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